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Abstract. We study sedimentation equilibrium of low-salt suspensions of binary 
mixtures of charged colloids, both by Monte Carlo simulations of an effective colloids- 
only system and by Poisson-Boltzmann theory of a colloid-ion mixture. We show 
that the theoretically predicted lifting and layering effect, which involves the entropy 
of the screening ions and a spontaneous macroscopic electric field [J. Zwanikken 
and R. van Roij, Europhys. Lett. 71, 480 (2005)], can also be understood on the 
basis of an effective colloid-only system with pairwise screened-Coulomb interactions. 
We consider, by theory and by simulation, both repelling like-charged colloids and 
attracting oppositely charged colloids, and we find a re-entrant lifting and layering 
phenomenon when the charge ratio of the colloids varies from large positive through 
zero to large negative values. 
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1. Introduction 

Suspensions of charged colloids are multi-component mixtures of mesoscopic (colloidal) 
particles and microscopic cations, anions, and solvent molecules. Due to the large size- 
and charge asymmetry between the colloids on the one hand and the microscopic species 
on the other, it is often difficult to treat all species on an equal footing in descriptions of 
e.g. the thermodynamic equilibrium properties of these systems. Given the asymmetry 
it is natural, and often practical, to view colloidal suspensions as " colloids-only" systems 
described by effective colloidal interactions in which the presence of the microscopic 
particles appear only through medium properties such as the dielectric constant e and 
the Debye screening length n~ l . For bulk systems in equilibrium one can in fact prove 
(by formally integrating out all the microscopic degrees of freedom in the partition 
function in a fixed colloid configuration) that such a course-grained description yields 
the exact thermodynamics; the problem is of course to actually perform the integration 
explicitly. Within linear screening theory it is possible to integrate out the ionic degrees 
of freedom explicitly, and the resulting effective colloid pair interactions are of the 
screened- Coulomb form oc exp(— nr)/r, with r the center-to-center distance between 
the colloidal pair £Q . Descriptions of suspensions of charged colloids on the basis of this 
effective pair potential have often been used over the years, see e.g. Refs.j3El for a few 
examples related to crystallisation of charged colloids, and Ref. 4J for a sedimentation 
study. 

However, in recent studies of sedimentation equilibrium of charged colloids at 
extremely low salinity this effective "colloids-only" view was not taken. Instead 
the experiments E], theory [TOJ [TTJ E2 EE31 d, and simulations [B] of 
sedimentation equilibrium were essentially all analysed within the so-called primitive 
model, in which the colloids and ions are treated as charged massive hard spheres and 
charged massless point particles, respectively, while the solvent is treated as a continuum 
with dielectric constant e. These studies have shown the existence of a spontaneously 
formed macroscopic electric field that lifts the colloids to much higher altitudes than 
to be expected on the basis of their mass. This electric field is the result of the 
competition between entropy (favouring a homogeneous fluid), electrostatics (favouring 
local neutrality), and gravitational energy (favouring the colloids at the bottom, but 
not the cations or anions because they are essentially massless) |13J. At high salt 
concentrations local neutrality can easily be fulfilled at low entropic cost because of the 
large number of ions, and hence the competition is the usual one between colloid entropy 
and gravity, giving rise to a barometric distribution at low packing fractions. At low 
salinity the situation is more complicated, since fulfilling local neutrality would imply, 
when the colloids are all close to the bottom because of gravity, that a large fraction of 
the counterions must also be close to the bottom. Such a state has a low ion entropy, 
because the ion distribution is far from homogeneous. It turns out, according to theories 
of e.g. Refs. jHllin], that the state with the lowest grand potential spontaneously sacrifices 
local neutrality at the bottom and at the meniscus, such that an electric field is generated 
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that lifts the colloids upward at the expense of gravitational energy, such that the colloids 
and the counterions can be rather homogeneously distributed over the available volume, 
thereby increasing their entropy. For colloids with buoyant mass m and charge Ze 
the generated electric field strength was shown to be E ±2 mg/Ze at low enough salt 
concentrations, such that the upward electric force ZeE and the downward gravitational 
force mg essentially cancel each other, leading to almost homogeneously distributed 
colloids. Here g is the gravitational acceleration, and e the proton charge. This lift 
mechanism involving a macroscopic electric field was confirmed by primitive model 
simulations of colloids and explicit ions ^3] , and by a direct measurement of the potential 
difference between the bottom and the meniscus [JJ. Moreover, quantitative agreement 
was found between the measured density profiles and the theoretically calculated ones 
in Ref. [Hj, while the measured profiles of Refs.jEl Ej agree with the theory of Ref. [2] 
which involves the electric field mechanism combined with charge regularisation. 

It is at first sight not at all obvious that these phenomena, which involve charge 
separation and macroscopic electric fields, can be described within the framework of 
the effective colloids-only picture. After all, in this picture the presence of anions and 
cations is only accounted for in the screening length k~ 1 that determines the range of the 
effective pairwise interactions. However, the theoretical study of Ref.[T3] showed that 
hydrostatic equilibrium between the external gravity field and the internal (osmotic) 
pressure of the effective colloids-only system (described by the simple and crude Donnan 
equation of state in Ref. [H|), actually predicts the existence of the electric field and the 
lift effect at low salinity, despite the underlying local density (and hence local neutrality) 
approximation. The predicted density profile is linear with height [T3] , not unlike those 
predicted by a density functional treatment of an effective colloids-only system with 
screened- Coulomb interactions jlj. Recently other accurate osmotic equations of state 
were employed and revealed the lift effect [To^ IT7] , and in fact Ref. shows that the 
colloid density profile that was experimentally measured and quantitatively described 
in terms of the primitive model in Ref. [S], can also be described quantitatively by a 
colloids-only theory with screened-Coulomb interactions. The inconsistency of a theory 
based on local charge neutrality that predicts an electric field without explaining its 
sources does, apparently, not lead to erroneous predictions for the density profile, at 
least not for monodisperse suspensions at the parameters considered. 

The question we address in this paper is whether the "colloids-only" picture can 
also account for the very recently predicted layering and segregation phenomena in 
equilibrium sediments of binary mixtures of charged colloids at low salinity [T%1 1191 12U] . 
These studies are actually direct extensions of the primitive model discussed before, but 
now with two (or more) colloidal components with different charges and masses. Not 
only was the lifting effect recovered in Refs. fTHllT^ I^Uj. but in addition segregation into 
layers was found, such that the two colloidal species order with height according to mass- 
per-charge: the colloids with the lowest mass-per-charge are found at higher altitudes. 
This implies the possibility of an inversion phenomenon, whereby sufficiently highly 
charged heavy colloids float on top of lighter colloids with a lower charge [THl H3 1201 • 
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In terms of the electric field this phenomenon can be understood, at least qualitatively, 
by first considering a pure system of colloidal species with mass mi and charge Zie, 
such that the electric field strength is E ^ m\gjZ\e according to (THj. Adding a trace 
amount of colloids with mass m 2 and charge Z 2 e does not change the field strength 
E, and hence the upward force Z 2 eE on the colloids of species 2 exceeds the downward 
force m 2 g provided m 2 /Z 2 < mi/Zi. In other words, if the mass-per-charge of the tracer 
species 2 is smaller than that of species 1, the former will float on top of the latter, even 
if m 2 > mx- I n the case that both colloidal species are present in substantial amounts, 
the segregation into layers was found with electric field strengths E^ ~ m^/Zje in the 
layer with species i [F£\. This implies a finite, nonzero charge density in the cross-over 
regime from one layer to another, as observed in Ref.|19j. The question we address in 
this paper is whether hydrostatic equilibrium (based on a local density and neutrality 
assumption) in a colloids-only system with screened- Coulomb interactions can catch this 
layering phenomenon. 

There is, in addition, a second motivation for our present study. It stems from 
recent experimental progress in preparing binary mixtures with oppositely charged 
colloids that do not aggregate but form equilibrium crystals [3j. It is of fundamental 
interest to investigate the consequence of opposite charges in the equilibrium sediment, 
in particular when there is a substantial mass difference between the two colloidal 
species: to what extent does the lighter species then lift the heavier one, and to what 
extent does segregation into layers take place. We therefore present not only results 
for colloidal mixtures with the same charge sign [Z\Z 2 > 0), but also with a different 
sign {Z\Z 2 < 0). Sedimentation in this latter parameter regime has, to the best of our 
knowledge, not been studied at all yet. We compare the results of computer simulations 
of the effective colloids-only system with theoretical calculations based on the Poisson- 
Boltzmann theory of Ref . ^H] ■ 

This paper is organised as follows. In section II we present the model of the 
suspension, and discuss it viewed as either a mixture of colloids and ions as in Ref.[T§] 
or as an effective colloids-only system that we simulate. In section III we present the 
resulting equilibrium density profiles, and section IV is devoted to a discussion. 

2. One model and two theories 

In this section we introduce the details and the notation of our model of a binary 
mixture of charged colloids that we use in this paper, and discuss two alternative ways 
of describing the equilibrium density profiles in sedimentation equilibrium. 

We denote the (fixed) colloidal charges, radii, and buoyant masses by Z^e, ai, 
and rrii for the two species i = 1,2, and assume that the charge is homogeneously 
distributed on the colloidal surface. The solvent is a continuum with dielectric constant 
e at temperature T, and the solvent volume is V — AH with A the (thermodynamically 
large) area in the horizontal plane and H the vertical height between the bottom of 
the sample (at height x = 0) and the meniscus (at height x = H). The gravitational 
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force points in the negative vertical direction, such that the potential energy of a colloid 
of species i at height x is given by niigx, with g the gravitational acceleration. It is 
convenient to introduce the so-called gravitational length Li = ksT jmig, which is the 
characteristic length scale of the barometric density distribution that holds in the dilute 
limit. 

The density profiles and the number of colloidal particles in the sample are denoted 
by Pi{x) and iVj, respectively, for species i — 1,2. We will also use the packing fraction 
profile rji(x) = (47r/3)afpi(x), and we characterise the density of the suspension by the 
overall packing fractions fji = (4ir/3)a^Ni/V = (1/H) J^ 1 dxrj^x). The average height 
of species i, or its center-of-mass, is defined as 



which can be seen as a rough indication of the nature of the colloidal profile. 

The suspension is imagined to be in osmotic contact with a reservoir of massless, 
monovalent cations and anions of charge ±e, with a total ion density 2p s , such that the 
Debye screening length (in the reservoir) is kT x = [Sti\bPs)~ 1 ^ 2 ■, where \b = e 2 jeksT 
is the Bjerrum length of the medium. Here ks denotes the Boltzmann constant. 

In the present calculations we restrict attention to spheres of equal radius, and we 
denote the (common) diameter of the spheres by 2a\ = 2a2 = o. 

2.1. Poisson-Boltzmann theory for the colloid-ion mixture 

It was shown in Ref.pH] that this system can be described by five coupled nonlinear 
equations for the following five unknown profiles: the colloid densities pi(x) and P2(x), 
the two ion densities p±(x) for the cations (+) and anions (-), and the electrostatic 
potential ip(x) or its dimensionless form <f)(x) = eip(x)/kBT. Two relations are given 
by the Boltzmann distributions p±(x) = p s exp[=p0(x)] for the ions, while the other 
three relations are a Boltzmann distribution for the colloidal densities and the Poisson 
equation for the potential, viz. 



where a prime denotes a derivative with respect to x, and where q represent 
normalisation constants such that A J Q dxpi(x) = iVj. We impose boundary conditions 
0'(O) = <p'{H) = 0, which corresponds to a globally charge neutral system without 
any external electric field. This set of five equations follows directly from the minimum 
condition on a grand potential functional that only involves ideal-gas contributions for all 
four species, a potential energy of the colloids due to their mass, and a Coulomb energy 
treated at a mean-field level ^H]- It is straightforward to solve this set numerically on an 
x-grid, for details see Ref-jTHj. From the solution for <f>(x) the magnitude of the electric 
field follows as (k B T/e)cf)'(x). 



J Q dxpi{x)x 
dxpi(x) 



(1) 




(2) 
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We note that these five equations are independent of the hard-core radii a, of 
the colloids, i.e. the short-range part of the (direct) correlations are not taken into 
account. As a consequence this level of approximation is not capable of describing 
packing effects properly. The radii dj are therefore only used to convert Pi(x) to rji(x). 
Despite this shortcoming the present theory agrees quantitatively with the primitive 
model simulations of Refs. [THJ |2T] . 

2.2. The colloids-only theory 

The Poisson-Boltzmann theory for the colloid-ion mixture differs, at least at first sight, 
substantially from a description of charged colloids based on the idea that colloids 
interact with each other through effective interactions that depend on the solvent 
properties such as the dielectric constant e, the Debye screening length and the 
temperature T. Standard linear screening theory P predicts that two colloids with 
charges Zie and Z^e and radii a, and aj, separated by a distance r, interact with the 
pair potential Vij(r) given by 



with the Boltzmann constant ks and the Bjerrum length Xb = e 2 /eksT as introduced 
above already. Here we ignore the attractive dispersion forces. Within this effective 
"colloids-only" description the presence of salt ions (with a concentration 2p s in the 
reservoir) is entirely included through the screening length We assume throughout 
that the screening length in the colloidal interactions is independent of the colloid 
concentration. Ignoring the counterion contribution to K will turn out to be entirely 
justified for all cases we study in this paper: the ratio y(x) = Y^=i ZiPi(x) /2p s is always 
such that \y(x)\ <C 1, i.e. the background electrolyte dominates the screening. Typically 
we find \y(x)\ ~ 0.01 — 0.2, and since the effective screening parameter can be written 
as k _1 (1 + y 2 )~ 1 ^ 4 according to the Donnan-like theory of e.g. Ref . |22[ l2"3*j . provided the 
packing fraction is not too high, we find with (1 + y 2 )^ 1 ^ — 1 — £/ 2 /4 that the effective 
screening length differs from the reservoir value k~ 1 by at most a percent. 

The effective "colloids-only" description based on Eq.(j3J) has proved to be very 
successful in describing many facets of colloid science, e.g. in the case of monodisperse 
suspensions it explains freezing into fee and bec crystals at sufficiently high densities 
[2], and it was recently used to describe crystal structures of oppositely charged colloids 
successfully 

In this paper we use Eq. (JSJ) in Monte Carlo simulations in a box of dimensions 
K x K x H, with K 2 = A the horizontal area. In all cases the vertical box height was 
taken to be H = 109<j, and the lateral width K = 10cr. We checked that K was large 
enough to exclude finite size effects. We employed periodic boundary conditions in the 
horizontal directions; in the vertical direction the system is bounded by hard walls that 
exclude the center of colloids at heights x < and x > H . 
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3. Results 

3.1. Heavy colloids 

We first consider binary mixtures of equal-sized colloids with a diameter a = 1950nm, 
a Bjerrum length Xb = 10.4nm, and screening constant such that kg = 1.2. This 
corresponds to a salt concentration p s = 2.4nM in the reservoir. In addition, we fix 
the colloidal charge of species 1 to Z\ — 76, and the gravitational length of species 1 
to L\ = 2.5a. These numerical values, which correspond to the experimental system of 
Ref. 0, give rise to an effective contact potential between the colloids of species 1 given 
by Vii(er) = 12/c^T. The mass of species 2 is taken to be larger than that of species 1, 
by a factor of 2, such that L 2 = 1.25a. 

In Fig^we plot the colloidal density profiles for three equimolar suspensions, with 
average packing fractions fji = fj 2 = 0.005, for three values of the charge of species 
2 given by (a) Z 2 = 300, (b) Z 2 = 0, and (3) Z 2 = -300. The box height is in 
all cases H = 109a. Case (a), for which Z 2 L 2 > Z\L\, shows a layer of the heavier 
species 2 floating on top of a layer of primarily species 1, both from the theory based 
on Eqs. (P?|) H5]. and from the colloids-only simulations using Eq.fjHJ). The agreement 
between the theory and the simulations is quantitatively reasonable though not perfect. 
In particular the packing effects shown by the simulations at length scales of the order 
of a few a close to the hard walls of the container are not reproduced by the theory, and 
the lifting effect for species 2 as predicted by the theory is quite a bit stronger than that 
found in the simulations. Nevertheless it is a striking observation that the colloids-only 
picture is capable of describing the layering phenomenon without invoking explicitly to 
a selfconsistent electric field. Case (b), where Z 2 = such that it represents a mixture 
of charged and uncharged colloids, shows a clear lifting effect for species 1, both from 
the theory and the simulations. The theoretical prediction for species 2 is a barometric 
distribution, which is bound to fail quantitatively in comparison with the simulations 
given the high packing fraction up to 0.5 at the bottom, which causes hard-sphere like 
density oscillations of species 1 close to the bottom in the simulations. Case (c) is for 
oppositely charged colloids, and shows again a lifting and layering effect, whereby a 
layer of the heavier species 2 floats onto a rather dense layer at x < 15a, both in the 
simulations and in the theory. Ignoring the density oscillations on a length scale of a 
few a in this bottom layer, the densities take the almost constant values ^i(x) ~ 0.035 
and r] 2 (x) ~ 0.018, i.e. there is no large systematic decay and the ratio of the two 
densities does not correspond to \Zi/Z 2 \. In fact, the agreement between theory and 
simulations in (c) is much less satisfactory than in the cases (a) and (b): not only are the 
packing effects in the bottom layer not captured by the theory (as before), but also the 
"smoothed" simulated profiles differ substantially from the theoretically predicted ones. 
We expect this to be caused by the strong correlations that build up between attracting 
pairs of oppositely charged colloids, with a contact potential of Vi 2 (a) = — 47.5k bT, 
and between the repelling pairs of species 2 with a contact potential V 22 (a) = 187/cbT. 
These correlations are not at all included in the theory. We checked explicitly that the 
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bottom phase in the simulations is fluid and not crystalline. 

An interesting observation that can be made on the basis of Figgis that species 2 
floats onto a relatively dense layer (with an almost constant composition) if Z 2 = ±300, 
while species 1 is floats onto a denser layer if Z 2 = 0. We can quantify this re-entrant 
layering effect by studying the average height hi as defined in Eq.((TJ) as a function of Z 2 , 
keeping Zi and all the other parameter fixed. The result is shown in FigJ21 and confirms 
the re-entrant phenomenon. It shows gross agreement between theory and simulation 
for all Z2/Z1, although the agreement is substantially better for like-charged colloids 
(Z 2 > 0) than for opposite colloids [Z 2 < 0). For Z 2 jZ\ > 2 the theory overestimates 
the lifting effect for the heavier particles and underestimates it for the lighter particles, 
in comparison with the simulations. For Z 2 jZ\ < —1 the theory predicts an essentially 
barometric distribution for the lighter species (with a lower absolute charge), whereas 
the simulations still show a considerable lift effect for species 1 in this regime, probably 
because they are correlated due to the attractions of the highly charged heavier colloids 
that float on top of them; this effect is poorly caught within the present theory. Note 
that for the equimolar case considered here the theory predicts h\ = h 2 for Z\L\ = Z 2 L 2 
(i.e. for Z 2 jZ\ = 2 here) and for Z\ ~ — Z 2 (since then <p{x) ~ 0). On the basis of 
the numerical results of FigJ2] these seem to be reasonable estimates for the crossover 
regimes. 

We also studied the effect of the screening length on the mean height of the colloids. 
The result is shown, for the parameters of Figda), in Fig|3| Here the theory is found to 
be not at all in agreement with the simulations, except perhaps for the decreasing mean 
height of the heavier species 2 with increasing kg. The crossover where h\ = h 2 is found 
to be at kg ~ 2.4 in the simulations, whereas the theory cannot locate a crossover at all 
for any reasonable kg. Moreover, the theory cannot catch the phenomenon of increasing 
hi with increasing kg for kg < 3. One could argue that the theory breaks down, not 
only quantitatively but even qualitatively, when kg > 2. 

3.2. Lighter colloids 

For the sake of comparison we repeated the study reported above to the case where both 
colloidal species are ten times lighter, such that L\ = 25a and L 2 = 12.5a. We left all 
the other parameters identical to those of FigsHEl 

In FigEJwe consider the three density profiles for the three charges (a) Z 2 = 300, 
(b) Z 2 = 0, and (c) Z 2 = —300, and find reasonably good agreement between theory 
and simulation for case (a) and (b), and now in fact also for (c). As in the case of the 
heavier colloids the theory still fails to describe the oscillating character of the profiles 
close to the bottom and the meniscus, but in all cases it can account for a "smoothed" 
version of the simulated profiles. We note that ) does not show a significant 

layering effect, in the sense that the density profiles of both species span the whole 
sample from bottom to meniscus. This can be understood from the theory of Ref. ^H] 
if one realises that both Z\L\, Z 2 L 2 ^> H, i.e. both species would prefer to be lifted to 
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much higher altitudes than the present system size H = 109cr allows. Unlike the layering 
effect, the lift effect remains clearly visible for this relatively small H. In Fig0fb) and 
(c) the layering effect does exist, although perhaps a bit weaker than in the case of the 
heavier colloids discussed above. We can again conclude that the lifting (and layering) 
effect can be accounted for within a colloids-only picture, and that there is a reentrant 
phenomenon with varying Z 2 from 300 through to -300 such that the heavier species 
is lifted provided that \Z 2 \ is large enough. This effect is quantified in FigJSJ which is 
the analogue of Fig|21 and shows the mean height hi as a function of Z 2 jZ\ at fixed 
Z\. The agreement between theory and simulation is more quantitative for these lighter 
colloids, at least for like-charged colloids. This is in line with our earlier notions that 
the differences between theory and simulation are probably due to correlations, which 
are weaker for lighter colloids since the system is then more homogeneous and hence 
more dilute at the bottom. The analogue of FigOl is shown in Fig® and shows the 
same poor agreement between theory and simulation as FigCJl the theory breaks down 
completely for na > 1.5, where it cannot even predict the correct ordering of the species 
with height. 

4. Conclusions 

We have studied the equilibrium sediment of binary mixtures of charged colloids in 
suspension, viewed both as a colloid-ion mixture within a Poisson-Boltzmann theory 
and as an effective colloids-only system with pairwise screened- Coulomb interactions 
in Monte Carlo simulations. The main conclusion is that the layering effect, whereby 
colloids order with height according to charge-per-mass as predicted by the Poisson- 
Boltzmann theory ^E1G2], can also be obtained within the effective colloids-only system. 
This result is not obvious, since the theoretical ion-colloid description invokes the ion- 
entropy and a selfconsistent electric field that pushes the colloids upwards against 
gravity, while the ions only occur very indirectly in the screening constant of the effective 
system. Nevertheless, these results suggest that both pictures are, at least qualitatively, 
different sides of the same coin. The inconsistency of the colloids-only system, where 
local electric neutrality is assumed whereas the (Donnan) potential varies with height 
and thus requires a non-vanishing local charge density at least somewhere in the system, 
is apparently no serious problem for describing the phenomenology quantitatively. 

It is of interest to try to understand the layering mechanism within the colloids- 
only picture qualitatively. For like-charged colloids one may argue that the colloids with 
the highest charge repel each other most strongly, such that they expand to relatively 
high altitudes, leaving the lower charged (and possibly lighter) colloids behind at lower 
altitudes. For oppositely charged colloids the situation is more subtle: the attraction 
between the most expanded highest charged colloids and the less expanded lower charged 
colloids is such that the tendency to layering is reduced. This effect depends, probably, 
sensitively on composition and charge, and has not been investigated here in any detail. 

At a quantitative level we can conclude that the theory does not catch the oscillatory 
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character of the density profiles that were found close to the bottom and meniscus in 
the simulations. This effect requires a better account of the short-range correlations, 
and will be addressed in future work. The theory is, all in all, in better agreement with 
the simulations for repelling like-charged colloids than for attracting oppositely charged 
colloids. This can probably also be traced back to the poor level of our theoretical 
treatment of the short-ranged correlations, which are more pronounced in the presence 
of attractions. Given the recent exciting new experimental developments in the study of 
oppositely charged colloids exhibiting equilibrium behaviour (as opposed to irreversible 
aggregation phenomena) j^j, there is good reason to attempt to improve the theory. 
Moreover, the theory is also shown to break down completely for kg > 1.5 — 2 (for the 
present choice of parameters), where short-ranged correlations become more important 
once more. Despite these shortcomings, our main conclusion should be that the layering 
phenomenon at low salt, say at kg J$ 1 — 2, is a real effect that can be understood from 
a colloid-ion mixture perspective as well as from a colloids-only perspective. We hope 
that these findings stimulate experimental study of these phenomena. 
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Figure 1. Simulated (full curves) and theoretically calculated (dashed curves) 
equilibrium sedimentation profiles for a binary mixture of lighter (species 1) and heavier 
(species 2) colloidal particles with charges Z\ = 76 and various charges Z 2 given by 
(a) Zi = 300, (b) Z2 = 0, and (c) Z2 — —300, as a function of altitude a: in a sample of 
height H = 109er. The common diameter of the colloids is a — 1950 nm, the mass ratio 
is m2/mi — 2, the gravitational lengths are L\ = ksT/mig = 2.5a and Li = 1.25<7, 
the total packing fractions are fji — 7)2 — 0.005, the Bjcrrum length is Xb = 10.4 nm, 
and the screening parameter is kg = 1.2. 
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Figure 2. Simulated (full curves) and theoretically predicted (dashed curves) mean 
heights hi of light (open circles, species 1) and heavy (solid squares, species 2) colloids 
as a function of the charge ratio Z2/Z1 for fixed Z\ = 76. All other parameters are 
identical to those of Fig^ 




Figure 3. Simulated (full curves) and theoretically predicted (dashed curves) mean 
heights hi of light (open circles, species 1) and heavy (solid squares, species 2) colloids 
as a function of the screening parameter kg for fixed Z\ — 76 and Z 2 = 300 . All other 
parameters are identical to those of Fig^ 
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Figure 4. Same as figd but for ten times lighter colloids with gravitational lengths 
L\ = 25a and L2 — 12.5a. 
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Figure 5. Same as fig|21 but for ten times lighter colloids with gravitational lengths 
L\ = 25a and L2 — 12.5<r. 




Figure 6. Same as figEl but for ten times lighter colloids with gravitational lengths 
L\ = 25a and L2 — 12.5a. 



